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ABSTRACT 

Rates of autoionization of H 2 near threshold is calculated within an internal 
conversion model. The present calculation remedies the apparent inability of 
the model to predict desirable rates as calculated by previous authors. Results 
are obtained which compare fairly well with the experimental estimates of 
Chupka and Berkowitz and alternative calculations of Berry and Nielsen for 
Av. f = 1 and 2. The results explicitly show the ^ dependence of rates for in- 
creasing principal quantum numbers n and follows the socalled 'propensity rule' 
for vibrational quanta involved. Unlike other theoretical calculations the present 
model provides a simple formula from which the rates may easily be calculated. 
Limitations and possible refinements of the model are discussed and generalisa- 
tions for diatomic molecules other than Hj, is indicated. 


AUTOIONIZATION OI- Uj NEAR THRESHOLD 
I. INTRODUCTION 

In molecular hydrogen if one of the electrons is removed sufficiently away from 
the nuclei, it may be considered to move in an approximate Rydberg State around 
the residual ion. High resolution experiments on photoionization and 
photoabsorption have not only confirmed the presence of such states for but 
also have yielded quantitative rates of autoionization from these states. The re- 
sults further reveal that the autoionization process, in fact, dominates over the 
photoionization process in the neighborhood of the threshold of ionization. Sev- 
eral theoretical calculations^'^ ® ’ have already been done to estimate the auto- 
ionization rates from these states. The calculations may be divided into two 
kinds of models. (Recently Ritchie® has calculated rates of autoionization within 
a semiclassical method, which are in qualitative agreement with experimental 
estimates for high vibrational states with single quantum transitions.) In the 
first kind of calculations, initiated by Berry^ the autoionization is assumed to be 
mediated by the nuclear kinetic energy terms which break down the Born- 
Oppenheimer separation by coupling the core vibrations with the electronic mo 
tion. We may call it the non-adiabatic model. In the second kind of model, first 
used for the present purpose by Russek et al.^ , the energy of vibrations is as- 
sumed to be mediated to the Rydberg electron directly through an internal con- 
version process, much used in nuclear problems before^ . 

In their first applications both the theories^- ^ failed to produce the estimated 
experimental rates, by falling short of by an order of magnitude or more. Sub- 

i 

sequently, however, the calculation of Berry and Nielsen^ using elaborate 
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non-adiabatic molecular calculations have improved the results to within a few 
times the experimental results. The experimental estimates^ may be uncertain 
within a factor of two or so and it may be considered that the non-adiabatic model 
is adequate enough to describe the autoionization process. However, due to the 
very nature of the theory, extensive numerical calculations are .lecessary to 
arrive at the estimates. The internal conversion iuodcl, on the ovher hand, is 
attractive for its simplicity and would be uso/ul if the model could be shown to 
yield the appropriate rates. 

From results of calculations wichin non-adiabatic models^ it has emerged that 
the monopole nuclear interaction is principally responsible for transmission of 
vibrational energy to the Rydberg electron. Previously it has been demonstrated 
by Russek et al.’ that there is an essential equivalence between the non-adiabatic 
and the internal conversion models. In view of these two facts we re-examine the 
internal conversion model and show in this paper that the model is capable of 
yielding the desired results via the so-called "0-0 transition" case well known in 
the analogous process of nuclear conversions^. 

n. MATHT.MATICAL FORMULATION 

Atomic units are used throughout the present formulation. We assume that the 
total Hamiltonian of the system can be iq>proximated by the sum of a Hamiltonian 
of the core and the Hamiltonian of the Rydberg electron, plus an effective inter- 
action between the Rydberg electron and the ion-core. Thus we write 

"... = «o(T) ♦ V.„(R,T) (1, 

where the core Hamiltonian may be given by 
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HcoreO^) 


^ E(R) 


( 2 ) 


1 ^ 1^2 

= _ V2 + — — 
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In (2) we neglect the rotational part of the Hamiltonian (we shall be 

I R2 

concerned solely with vibrational autcionizatiors for transitions between ground 
J = 0 rotational states). 


In the present first order theory we further confine ourselves to the ground 
electronic state of the core and neglect all other effects of the core electron 
besides its effect in generating the groimd state potential energy £ (K) (which 
supports the nuclear vibrations). 


The Hamiltonian of the Rydberg electron is defined by 






( 3 ) 


Thus the autoionizing states are characterized by the usual hydrogenic wave 
functions with quantum numbers nlm. 


We, therefore, have the potential interaction 


V(R,r) 



( 4 ) 


where the first two terms represent the interaction of the Rydberg electron witii 
the nuclei (with charges and Zj) and the third term is the correlation inter- 
action with the core-electron at F . The last term - should be included in 

r 

equation (4) in order that the definition (3) of the Rydberg states is preserved 
throughout. It may be noted in this connection that in any decomposition of the 
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total Hamiltonian, a part of the Hamiltonian that Is included in defining the 
'unperturbed' states of the system must not be present in the interaction (part) 
that mediates transitions between these states. By expanding (4) in Legendre 
Polynomials it is easily seer, that for large r, there is no monopole term in the 
expansion and that the coupling between the nuclei and the Rydberg electron is 
of the quadrupole order (the dipole interaction disappears due to the fact that 
the core is homo-nuclear in nature). It is this quadrupole interaction 

(R/2)2 P,(R • f) 

which was considered by Russek et al.^ to be responsible for autolonizing 
transitions. The contribution from the quadrupole interaction, however, turned 
out to be too small by one or two orders of magnitude^. One notes now that 
there is a second part of the interaction (4) that corresponds to r < R/2. Thus 
neglecting the correlation term (as assumed in the present model) and expanding 
we find 


V(ll,r) 



when r < R/2. 


(5) 


We shall consider only the leading terms in (S). Although at first sight it ap- 
pears that the nuclear monopole terms being independent of coordinates of the 
Rydberg electron is incapable of producing the coupling required to mediate the 
transition: on closer examiiuUions however, one finds the situation to be different. 
Thus we shall find that to the extent the 'tail* of the Rydberg electron penetrates 
the core and overliqps with the nuclear wave function, the monopole interaction 
provides a non-vanishing contribution to the transition matrix elements. 


According to our decomposition of the total Hamiltonian we write down the first 
order wave function of the whole system to be 

'/'i(K.r) = X,, (R) (r) Yp^(0 (6) 

where (R) are vibrational wave functions of the core, R„p (r) is the radial 
wave function of the Rydberg electron, and Y£ ^ (r) are the normalized spherical 
harmonics for the angular wave functions of the electron. (In the present treat- 
ment we are interested in the vibrational autoionization process itself and 
therefore, have explicitly avoided introducing rotational wave functions of the 
core. In future we intend to publish the corresponding theory for rotational 
autoionizations involving ^propriate frame-transformation properties of the 
system. At least one case of rotational autoionization in H 2 has already been 
observed^). In the final state, likewise, we assume the wave function to be 

'/'fCR.n = X,,(R)<Pp(r) (7) 

where (?) is the coulomb wave for the ionized electron moving in the direction 
of the final s^ate momentum tc:*° 

<Pf(T) = ^ e R|,f (r) Y^i^, (r) Y£*|j^, (k); agi = coulomb phase. (8) 

f-o 

The functions <P|- (r) are normalized in momentum space by the relation 

d7 = S(k-k') (9) 

The transition probability for autoionization may now be given by*: 
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P, ■> \ ^ C(W) 


( 10 ) 



where l{k) is the density of final states, which, in atomic units, is equal to the 
momentum k . The above ejq;)resBlon corresponds to the probability of ejection 
in the differential solid angle dO^. Substituting the wave functions (6) and (7) 
and the interaction (5) in equation (10) and performing the angular integrations 
over dr, we obtain 


dT 

dn._ 


= 2^ jjdK X,, (R) (-|)x,,(R) X 'r f(r) R,^(r) r^dr 


( 11 ) 


j*dRX,,(R) (7)^f 


(r) r* dr 


(iff.r; v,-.(k)) 


We note that the penetration integrals over r in (11) range within a small region, 
between 0 to R/2. Thus the radial functions may be conveniently iq>proximated 
by‘°: 


R„f(0 


_ 1_ / (n +f )! I'* , 

■ (2f + 1)! * \(n-^ - 1)!/ " 


( 12 ) 


Since all electrons considered are ejected with energy < < 1 a.u., we may also 
iq)proximate R,,f (r) by the leading term of its small energy «q)ansion^: 


Lim^.o (*■) % ♦ I (/i7) 


(13) 
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Substituting (12) and (Id) in (11) and integrating over r , we find: 




(n > f ) ! ^ 1 

(n - f - 1) {(2P + 1) !}^ 




The integrations over r are performed analytically by noting that * ' 


I dz z*' = z*' J^,(z) . 


The total autoionization probability is readily obtained by integrating (14) over 
all angles of ejection: 

T = 277 • A„P • (a.u. )■> (17; 

To obtain the rate per second we need only to divide (15) by the a.u. of time^^, 
to = 2.4189 X 10" sec. 


in. RESULTS AND DISCUSSIONS 


From e}q)res8ion (15) and (17) we readily find that the rates of autoionization 
decrease as — ; with increasing n (as was noticed by Bardsley*'*|. 

The vibrational integrals | v, *• ^ | ^ (which are independent of n ) are evaluated 
numerically, using Morse functions (see Appendix) for the ion-core. 
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In Table I we compare our calculated rates with the theoretical results of Berry 
and Nielsen^ and tlie experimental estimates of Chupka and Berkowitz^. We 
also quote for comparison purposes the contribution from the quadrupole term 
as calculated by Russek et al.’ 

In Table II we present some results corresponding to vibrational transitions 
involving two quantum transitions = 2, and compare them with the available 
theoretical results of berry and Nielsen^. We find in agreement with these 
authors that the rates for Ai/. , = 2 are much smaller than those for Ai/j, = 1, 
thus satisfying the 'piopensity rule* pressed by Berry 

In Table in we present some results for high vibrational transitions with , = 

3 and 4. We find that unlike the At/ . ^ = 1 and 2 cases these results differ 
greatly for most cases with the results of Reference 7. = 4 v, =0 case 

seems to be an exception I) Although these rates continue to satisfy the 'pro- 
pensity rule' the disagreement may well be due to sensitive dependence ol multi- 
ple quantum transitions on the exactness of vibrational wave functions (which we 
have iq^roximated here simply by Morse functions). 

We note that formula (17) does not distioguish between a and rr states, but can 
yield rates of autoionUation from states other than p-states ^ populated by 
non-optical means). We also have noted that if one omits the term j from the 
^fective interaction (5) the rates turn out to be somewhat higher (as may be 
expected) and tend to agree more closely with the experimental rates given in 
Table 1. However (as explained in Section II), for a satisfaotmry deftoitlon cf 
the Rydberg states (equation (3) ) we have retained the ^ term in the effective 
interaction (equatimis (4) or (5) ). 
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In a more detailed treatment (which may be necessary if the experimental 
determination of rates improve significantly in fugure) one could also take ac- 
count of the core-electronic functions and the electron-electron correlation 
explicitly. We expect this refinement to reduce (due to the repulsive nature of 
electron-electron interaction) the present estimates to a certain extent, but 
leave the essential feature of the model unaffected. 

Within the approximations introduced in the present paper the generalization of 
the method to diatomic molecules other than Hj (e.g., Nj, CO, etc.) would con- 
sist simply in substituting appropriate vibrational functions (or Morse functions), 
for the respective molecular iont, in expression ( 16 ) for the vibrational matrix 
elements. 
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APPENDIX 


The vibrational wave functions, x:,,(R), for diatomic molecules are given (usually 
quite accurately) by Morse functions*^- 

X,(R) = N, e‘«« F^(2) 

where 
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Table I 


Autoionization rate (see* ^ ) of np\ Rydberg state v,f Hj. Av. ^ = 1; CB; Ex- 
perimental estimates of Chupka and Beckowitz^, BN; Theoretical calculations 
of Berry and Nielsen' and RPB: Theoretical calculations of Russek et al.^ 
(quadrupole contribution), PC: Present Calculation (the rates forcr and n states 
are not distinguishable in the present model) 






3 2 



nf k 

CB 

PC 

Sp.T 

4.0 X 10 

1,9 X 10** 

8P7T 

(3.0 X 1C‘^) 

1.9 X 10** 

9p7T 

(3.0 X 10**) 

1.0 X 10** 

8po- 

1.3 X 10*2 

5.0 X 10** 

8p7r 

(9,0 X 10**) 

1 

5.0 X 10** 

b 

a 

00 

2,7 X 10*2 

9.6 X 10** 

ISpcr 

(1.3 X 10*2) 

1.2 X 10** 

7 per 

(9.0 X 10**) 

2.3 X 10*2 

(8 per) 

(5.0 X 10*2) 

1.6 X 10*2 

9pCT 

1 

2,7x10*2 

1.1 X 10*2 

1 

lOper 

(3,0 X 10*2) 

8.4 X 10** 

13pa 

(3.3 X 10*2) 

3,8 X 10** 

8p 

2,0 X 10*2 

2,6 X 10*2 

9pa 

(lOper) 

5,8 X 10*2 
(1.7 X 10*2) 

1,8 X 10*2 
1,3 X 10*2 











1,6 X 10*2 

1,02 X 10*0 

1.1 X 10*2 



8.0 X 10** 

0.53 X 10*0 
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